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ABSTRACT 

We investigate the lifting effect of rotation on the thermal evolution of 
carbon-oxygen white dwarfs accreting carbon-oxygen-rich matter. We find that 
rotation induces the coohng of the accreting structure so that the delivered 
gravitational energy causes a greater expansion with respect to the standard 
non-rotating case. The increase in the surface radius produces a decrease in the 
surface value of the critical angular velocity and, therefore, the accreting white 
dwarf becomes gravitationally unbound (Roche instability). This occurrence is 
due to an increase in the total angular momentum of the accreting white dwarf 
and depends critically on the amount of specific angular momentum deposited 
by the accreted matter. If the specific angular momentum of the accreted 
matter is equal to that of the outer layers of the accreting structure, the Roche 
instability occurs well before the accreting white dwarf can attain the physical 
conditions for carbon burning. If the values of both initial angular velocity 
and accretion rate are small, we find that the accreting white dwarf undergoes 
a secular instability when its total mass approaches 1.4 Mq . At this stage, 
the ratio between the rotational energy and the gravitational binding energy 
of the white dwarf becomes of the order of 0.1, so that the star must deform 
by adopting an elliptical shape. In this case, since the angular velocity of the 
white dwarf is as large as ~ 1 rad s^^, the anisotropic mass distribution induces 
the loss of rotational energy and angular momentum via gravitational wave 
radiation. We find that, independent of the braking efficiency, the white dwarf 
contracts and achieves the physical conditions suitable for explosive carbon 
burning at the center so that a type la supernova event is produced. 

Subject headings: stars: accretion, rotation - supernovae: general - white dwarf 
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1. Introduction 

Type la Supernovae (SNe la) are one of the most powerful tools for determining the 
rate of expansion and elucidating the geometrical structure of the Universe. They are 
intrinsically very luminous objects and, provided the empirical relation between the absolute 
magnitude at maximum light and the rate of decline from maximum (Phillips 1993) proves 
to be theoretically well founded and universally valid, they may eventually be used as fairly 
reliable distance indicators. For each observed event, it may then be possible to estimate 
reliably the luminosity at maximum and, hence, to determine the distance modulus. 

Despite the pivotal role of the observed correlation in observational cosmology, a 
theoretical justification for the correlation does not yet exist. Nevertheless, a Phillips 
relation which has been constructed from a sample of only nine local {z < 0.5) supernovae 
(SNe) has been used to estimate the distances to many more SNe at much larger z, assuming 
that no differences at all exist between events produced by different stellar populations in 
the past. 

Finally, from the theoretical point of view, a definitive identification of the progenitor 
system and an understanding the full nature of the explosion mechanism are still missing. 
At present it is accepted that SNe la are produced by the thermonuclear disruption of 
a carbon-oxygen (CO) white dwarf (WD) which has attained the Chandrasekhar mass 
{Mch ). Hoyle & Fowler (1960) were the first to examine the nucleosynthetic characteristics 
of such an explosion, without specifying how the initial configuration comes about. Arnett 
(1969) demonstrated that, without the mass loss which in real hfe terminates its evolution, 
an intermediate mass single star would ultimately develop a CO core of Chandrasekhar 
mass. We now know that single intermediate mass stars end their nuclear burning evolution 
with the ejection of their hydrogen-rich envelopes (illuminated for a brief time as planetary 
nebulae) and evolve into white dwarfs of mass considerably smaller than the Chandrasekhar 
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mass. The consensus view has become that the progenitors of type la supernovae are 
close binary star systems. The most traditional scenario is the single degenerate (SD) 
one in which it is assumed that the accreted matter is supplied by a companion which is 
composed of hydrogen-rich matter (Whelan & Iben 1973). An alternate scenario is the 
double-degenerate (DD) one which argues that evolutionary paths exist for producing 
two CO white dwarfs which are close enough to merge in a Hubble time (due to angular 
momentum loss by gravitational wave radiation) and whose combined mass exceeds Mch ■ 
This second scenario, examined by Iben & Tutukov (1984) and by Webbink (1984), has been 
questioned from time to time. Two main features make the DD scenario attractive: 1) the 
total absence of hydrogen in the progenitor system and 2) the fact that the more massive 
WD approaches the critical mass in a very natural way. Yet, in spite of huge observational 
efforts, only one good candidate SN la progenitor has been identified (KPD 0422-1-5421 - see 
Koen, Orosz & Wade 1998); in fact, this is the only system yet found the total mass of which 
is of the order of Mch and the components of which will merge in less than a Hubble time 
due to gravitational wave radiation (GWR). Other observed systems with an orbital period 
short enough for a merger to occur in less than a Hubble time have a total mass which is too 
small for the system to qualify as a SN la progenitor (see Saffer, Livio & Yungelson (1998) 
and also Livio (2000) for a recent review). Nevertheless, population synthesis estimates 
suggest that enough massive close binary systems are produced in nature to account for the 
observed frequency of SNe la events (Iben, Tutukov, and Yungelson (1997); Livio (2000), 
and references therein) . Most of these systems are undetectable because of their very low 
intrinsic luminosity. However, thanks to an increase in luminosity due to viscous heating 
(Webbink and Iben 1988; Iben, Tutukov, and Fedorova 1998), perhaps a very few can be 
seen just before merger, due to an occultation-modulated light curve. 

It is commonly understood (e.g., Tutukov and Yungelson (1979), Iben and Tutukov 
(1984), Benz et al. (1990), Rasio & Shapiro (1995), Mochkovitch, Guerrero & Segretain 



- 5 - 



(1997)) that if two CO WDs merge, the less massive component disrupts into an inflated 
disk around its companion. In this situation it has been argued that, due to a high effective 
accretion rate from the disk, an off-center carbon ignition occurs and the more massive 
star will be converted into an oxygen-neon-magnesium (ONeMg) WD (Nomoto & Iben 
(1985), Saio & Nomoto (1985), Kawai, Saio & Nomoto (1987), Timmes, Woosley & Taam 
(1994), Saio & Nomoto (1998)). These discussions do not take into account many of the 
complexities involved in the merger process (see, e.g., Iben (1988)). Perhaps one of the 
largest uncertainties plaguing the one dimensional simulations which attempt to address 
the problem is the unknown rate at which angular momentum diffuses out of the thick 
disk, allowing matter from the disk to settle onto the white dwarf and adopt a structure 
describable by a traditional one-dimensional model. 

In this paper, we first present results of numerical experiments in which CO-rich matter 
is accreted radially by an initially "cooled-down" CO WD of 0.8 Mq at rates in the range 
10~^ — 10~^Mq yr~^. We explore in detail the thermal response of the accreting structure 
and point out the main physical parameters which drive the accretion process and determine 
the final outcome. These experiments provide a reference point for understanding results of 
further computations in which the lifting effect of rotation is taken into account. In choosing 
an initially cooled-down structure, we ignore the possibility that viscous heating prior to 
merger may be significant, but this is legitimate since there is at present no theoretical or 
experimental framework for estimating the magnitude of the relevant turbulent viscosity. 

The plan of the work is as follows: in §2 we present the input physics and the adopted 
numerical techniques; in §3 we discuss the evolution of the accreting structure as a function 
of the accretion rate for standard (non-rotating) models; in §4 the hfting effect of rotation 
on the evolution of accreting models is discussed. Conclusions are summarized in §5. 
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2. Input Physics 

The accretion experiments have been computed with an updated version of the 
FRANEC code (Chieffi & Straniero 1989), the main differences being discussed by Cassisi, 
Iben & Tornambe (1998). The equation of state is that computed by Straniero (1988) with 
succeeding upgrades; for the radiative low temperature opacities (T < 12000 K) we adopt 
the tables provided by Alexander & Fergusson (1994) whereas for higher temperatures the 
tables provided by Huebner et al. (1977) are used. The contribution of conductive opacities 
has been included following the prescriptions by Mitake, Ichimaru & Itoh (1984) and Itoh, 
Hayashi, & Kohyama (1993). The accretion process has been computed assuming that the 
accreted matter has the same specific entropy as the matter at the surface of the WD (for 
a discussion of this topic see Taam (1980a,b), Iben (1982), Nomoto & Iben (1985)). We 
have adopted values of the accretion rate in the range 10~^ ^ 10^^ M© yr~^; the upper limit 
is an estimate of the maximum accretion rate at which the transferred matter is accepted 
by the accreting WD (Eddington limit), whereas the lower limit has been fixed so that the 
accretion process can produce a SN la event well before the onset of crystallization in the 
inner zones of the WD. The chemical composition of the accreted matter has been chosen 
to be equal to the chemical composition of the matter in the surface layers of the CO WD, 
namely Xi2c = 0.507, Xwq = 0.482 and X20Ne = 0.011. 

EDITOR: PLACE FIGURE |1| HERE. 



3. Accretion with Standard Conditions 
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3.1. The Initial Model 

The initial model is a 0.8 Mq CO WD obtained by evolving a solar metallicity pure 
He star from the helium main sequence up to the cooling sequence (for more detail see 
Limongi & Tornambe (1991)). In the computation of the central helium-burning stage, 
semiconvection has been accounted for but breathing pulses at the border of the convective 
helium core have been neglected, as suggested by the results of Caputo et al. (1989). 

When the model reaches the cooling sequence and the helium-burning shell switches off, 
the oxygen profile in the core exhibits a maximum and the carbon profile exhibits a minimum 
(see panel a in Fig. 1). This particular feature forms when, after the exhaustion of helium 
at the center, the radiative helium-burning shell goes through the region where helium was 
partially depleted during the outward extension of the convective helium core in the central 
helium-burning stage; in this zone, due to the reduction in the local helium abundance, 
a-captures on ^^C compete with 3a reactions. The extension in mass of this oxygen-enriched 
region depends on the efficiency of semiconvection, while the level of the peak depends 
on the adopted cross section for the ^^C(q;,7)^^0 reaction. Adopting the prescriptions in 
Salaris et al. (1997), the chemical profiles have been modified as in panel b of Fig. 1 and 
then the cooling phase has been computed. The physical characteristics of the cool model 
are: log(L/L0) = -2.858, log(reff) = 4.049, log(pc) = 7.036, and log(Tc) = 7.015. 



3.2. Time Scales of Relevance for an Accreting White Dwarf 

The thermal evolution of a CO WD accreting CO-rich matter can be discussed in terms 
of several relevant time scales. For example, the final outcome of the accretion process 
can be understood by examining the behavior in time of the ratio of a thermal diffusion 
time scale (Henyey & L'Ecuyer 1969) and a compressional heating time scale (Sugimoto & 
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Nomoto 1977, Nomoto 1982a) 

Ignoring all but gravothermal energy generation, conservation of energy requires that 

On the left hand side of this equation, U is the thermal energy per unit mass, P is the 
pressure, V is the specific volume {V = 1/p), p is the density, and derivatives with respect 
to time are co-moving derivatives in Lagrangian coordinates. On the right hand side, 
assuming that energy flow is by radiation and conduction only, 

IGaT^dT^ 

where T is the temperature, a is the Stefan-Boltzmann constant, k is the opacity, and f is 
a unit vector in the radial direction. 

A local compressional heating time scale can be defined by setting the left hand side of 
eq. (1) equal to zero, giving 

1 /I dU\ f P \ d\np 



Tcomp \U dt J ^.^^^^.^ \pUj dt 



(3) 



Since the quantity ( — ] is nearly constant, varying in value from 1 (electrons with 

2 

non-relativistic velocities) to - (relativistically degenerate electrons), it is reasonable to 

,1.1 .Id Inp 
adopt the simpler expression = — - — . 

^comp "'^ 

Nomoto & Sugimoto (1977) set 

1 d Inp / d lnp\ / d Inp \ d\nM 



(d lnp\ _ / (9 Inp \ d\nM 
\ dt ) ~ \dhiMrj^ dt ' 



^COmp dt y J ^ y L^iiiii-j.^ y ^ 



(4) 



where M is the total mass, M^. is the mass at radius r, and differentiation is with respect 
to mass fraction q = and time t as independent variables. Neglecting the first term in 
eq. (4), Nomoto & Sugimoto adopt 

1 _f d\np\ d InM _ GM^ ( d \np\ 1 

^ ~ \dhiMr]^ dt ~ Anr^P [d InpJ^ ^ ^ 
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where 

(d\YiM\-^ 

(6) 



dt 



is the time scale for accretion. 



Differentiating with respect to q and M as independent variables, Nomoto (1982a) 
gives 



d Inp 



/ d Inp \ I d lnp\ 
showing that, in all approximations, Tcomp is proportional to Tacc- 



A local thermal diffusion time scale can be defined as the time required by a thermal 
signal to propagate over a specified distance. Equations (1) and (2) give 

(dU r,<^^\ _ ^ dT (d\np\ dP_(d 2\(l6aT^dT\ 

dP 2 

where Cp is the specific heat at constant pressure. If we set —— — 0, neglect - relative to 

3 ^ 

— , and assume that — varies much less rapidly with distance than — , then 
dr Kp dr 

dT d^T 



where 



, 16cT , , 



Equation (9) is the well known equation for the diffusion of heat and its solution for 

many situations is proportional to (e.g., Margenau & Murphy 1943) 

where 5r is the distance from an initial delta-function perturbation in T . This motivates 
the choice of 

Tdiff,5r- - 4^ (12) 
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as a measure of, or "time scale" for, the diffusion of heat across the distance 5r. 

As suggested by Sugimoto (1970), an appropriate length scale is the local pressure 
scale height Hp, leading to a "local" diffusion time scale 

Hp 

Tdiflf.loc = (13) 

(Sugimoto chooses a value 4 times larger than this). One might naively assume that a time 
scale for the propagation of heat over many scale heights is given by the sum of local time 

scales over these scale heights. However, as discussed by Henyey & L'Ecuyer (1969), the 
appropriate approximation is instead 

fri 



(14) 



raiff^iob = ^dr) = — j^^ (^) pdr 
where the integral extends from ri to r2- 



In the following, we will make use of a "global" compressional heating timescale defined 
as Tch = Tacc s-nd of a "global" diffusion time scale T^is equal to the value of Tdifr,giob given 
by eq. (14) when ri — and r2 = -R, where R is the radius of the star. 

The gravitational energy released by accreted matter is partly stored in layers near 
the surface as thermal energy, but partly diffuses inward, increasing the thermal content 
of layers ever deeper below the surface. If Tdifr > Tch, heating of layers near the surface 
may lead to the off-center ignition of carbon, well before the white dwarf reaches the 
Chandrasekhar hmit (Nomoto 1982a,b, Nomoto & Iben 1985). If Tdifr < Tch, energy can be 
transferred inward efficiently enough that carbon ignition occurs at the center or, if oxygen 
has replaced carbon at the center, at the base of a carbon-rich layer (Isern et al. 1983, 
Hernanz et al. 1988). 

Choosing D = R^u, ~ lO'^i?©, T = Tc ~ lO'^ K (isothermal structure), Cp ~ 10^ erg/g 
K, p equal to the average density of the structure (~ 10^ g/cm^), and k equal to the mean 
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conductive opacity (/t ~ 10~^ cm^/g), our initial model is characterized by rjifj ~ 3 x 10'' 
yr. Insisting that T^i^ = r^h, the maximum value for the accretion rate suitable for the 
efficient transfer toward the inner zones of gravitational energy released by accretion can be 
estimated as Mmax ~ 3 x 10^^ Mq yr"^. It follows that, if M < Mmax, the thermal energy 
excess in the outer layers is transferred to the interior and carbon burning is ignited at the 
center. 

The value of Tdifr just estimated refers to the initial model. Since Tdifr actually evolves in 
time as a consequence of the evolution of the physical conditions in the accreting structure, 
it is necessary to examine its value in the evolving models. 

EDITOR: PLACE FIGURE | HERE. 

In Fig. 2, the evolution of T^m (solid lines) is shown in comparison with the evolution 
of Teh (dashed lines) for all of the models computed in the present work, as listed in Table 
1. It is evident that, for M < M^ax {panel i), the mass of the accreting WD reaches the 
Chandrasekhar limit, so that central carbon burning occurs. It is apparent also that, for the 
cases with Mmax < M < 10^^ Mq jt~^ {panels etoh), the gravitational energy deposited 
in outer layers is transferred efficiently toward the center since, during the accretion process, 
the compressional heating time scale increases while the thermal diffusion time scale 
decreases until the ratio of the two time scales is reversed. This occurrence can be easily 
explained by observing that, on increasing the total mass, the radius and the conductive 
opacity of the accreting structure decrease while density and temperature increase. This 
behavior is common to the entire set of computed models, but if M > 10~^ Mq yr~^ 
, carbon-burning ignition occurs in the outer layers well before a significant amount of 
the deposited gravitational energy can be transferred toward the center. Fig. 2 discloses 
further that, for small enough accretion rates (say, less than 5 x 10"'' Mq yr~^), the thermal 
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diffusion time scale profile has two different slopes, the inflection point corresponding to the 
time when the thermal wave reaches the center. 

EDITOR: PLACE FIGURE | HERE. 

This behavior can be better appreciated in Fig. 3 where the evolution of the mass 
coordinates of the maximum temperature (M-r-max) and of the central temperature (Tc) are 
reported together with rjifj as functions of the total mass of the accreting structure for the 
case M = 5 X 10^^ Mq yr^^. Where the central temperature increases most rapidly with 
mass (mass coordinate ~ 1.05 ± 0.05 Mq), the rate of change in the slope of Tdifr decreases 
rapidly, becoming almost fiat. The maximum in temperature does not reach the center 
until the whole structure has become isothermal. Thereafter, the evolution of the central 
temperature is driven by homologous compression. 

The final sudden increase of the central temperature and the corresponding decrease in 
TdifT when the total mass of the accreting WD reaches ~ 1.35 Mq is a consequence of the 
strong compression due to the approach to the Chandrasekhar mass. 

A discussion of global behavior can also be conducted in terms of local time scales. 
By way of illustration, we examine the cases of accretion at rates marking the boundaries 
of the region we have explored, namely, 10~^ Mq yr^^ and 10^^ Mq yr~^. In Fig. 4, the 
evolution of the temperature, the local compressional heating time scale rch,ioc given by eq. 
(2), and of two ratios involving Tch^ioc , Tdis.ioc given by eq. (11), and Tevo, the time elapsed 
since the beginning of the experiment, are shown for the case 10~^ Mq yr~^. 

EDITOR: PLACE FIGURE | HERE. 

In the inner regions (M < 0.65 Mq), the heat diffusion time scale is long compared 
with the compressional heating time scale, so the increase in temperature is due entirely to 
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local compressional heating. In central regions (0.65 < M < 0.8 M©), the heat diffusion 
time scale decreases in the outward direction to become comparable to the compressional 
heating time scale, and heating is due to a combination of compressional heating and inward 
diffusion, the contribution of diffusion increasing outward with respect to the contribution of 
compression. In the outermost accreted layers (M > 0.8 Mq), Tdiff^ioc is everywhere smaller 
than Tchjoc and, since also everywhere Tdifr^ioc < Tcvo, diffusion efficiently redistributes energy 
liberated by compression both inward and outward. The peaks in the temperature profiles 
show that the "watershed" marking the place where the heat ffow changes directions is 
approximately at the center of the accreted layer. 

Figure 5 shows the evolution of the same quantities displayed in Fig. 4, but for the case 
10^^ Mq yi^^. In this case, Tdig^ioc is always larger than Xchjoc everywhere in the interior; as 
Tevo increases, compressional energy released everywhere is redistributed in such a way that 
the whole structure becomes nearly isothermal. 

It is worth pointing out that, over the matter in the initial white dwarf, the mass 
averaged value of Teh is essentially race, whereas, in the accreted layer, Teh is everywhere 
much less than Tacc- So, what we really do when we compare rdig^giob with race is compare 
the time it takes energy deposited in accreted layers to make its way into the deep interior 
of the original white dwarf with the time it would take these layers to heat up without 
inward diffusion. 



EDITOR: PLACE FIGURE | HERE. 
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3.3. The Accretion Experiments 

In Fig. 6 we plot the evolution of the accreting structures in the HR diagram while 
in Fig. 7 we report the evolution of the surface radius as a function of the total mass 
for all the models listed in Table 1. It is evident that, at the beginning of the accretion 
process, the structure returns along the cooling sequence discussed in §3.1 and expands 
since the gravitational energy deposited at the surface causes a local increase in temperature 
and, hence, in luminosity. It is important to keep in mind that the initial WD is a very 
compressed isothermal structure so that the local efficiency of the inward thermal diffusion 
is very small [rdm is very large); indeed, all of the energy delivered by the accreted matter 
is locally stored as thermal energy. The increase in pressure related to the increase in 
temperature leads to a rapid expansion of the outer layers; the expansion induces a decrease 
in the local density while the continuous deposition of matter produces an increase in 
the local temperature. In fact, the value of T^is in surface layers begins to decrease and, 
when the energy excess can be transferred to the underlying zones, that is, when the local 
Tdifj becomes less than Teh, the structure begins to contract. As is clearly shown in Figs. 
6 and 7, the larger the accretion rate, the larger is the radius at the end of the expansion 
phase and the larger are the luminosity and the effective temperature; as the accretion rate 
is increased, the rate at which gravitational energy is deposited in surface layers is increased 
and, hence, the induced expansion becomes larger. 

EDITOR: PLACE FIGURE ?? HERE. 

EDITOR: PLACE FIGURE HERE. 

In Fig. 8, we plot the evolution of the center of the accreting WD in the p — T plane 
and in Table 1 we report for each adopted value of the accretion rate the final mass (in 
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Mq), the mass coordinate (in Mq) at which carbon burning is ignited, and the value of 
temperature (in K) and density (in g cm~^) at this point. Only for M smaller than 10"^ 
Mq yr~^ does central carbon ignition occur. Three different phases can be distinguished in 
the thermal evolution of models which ignite carbon burning at the center: 

EDITOR: PLACE FIGURE | HERE. 

• the first phase, during which the thermal evolution of the center is determined by 
homologous compression. In this phase, the increase in mass of the WD induces the 
compression of the whole structure while the gravitational energy delivered by the 
accreted matter is stored in the outermost layers; 

• the second phase, characterized by a rapid increase in the central temperature. This 
rapid increase in central temperature is a consequence of the fact that the gravitational 
energy initially localized in outer layers has diffused to the center. During this phase, 
homologous compression continues to be active and drives the evolution of the central 
density. This phase is concluded when the structure becomes almost isothermal; 

• the third phase, during which thermal diffusion and homologous compression are 
both effective. At this stage, the structure is isothermal, as it was when the accretion 
process began, but now the temperature level is higher; as a consequence, gravitational 
energy is not stored for long in outer layers but is almost immediately transferred to 
underlying zones (see the drop in T^is in Fig. 3 from the first to the third phases). In 
addition, the gravitational energy delivered by the accreted matter is negligible with 
respect to the heating produced by homologous compression, the efficiency of which 
is independent of the thermal content of the accreting structure; as a consequence, 
during this phase, homologous compression becomes the leading factor which governs 
the increase in the interior temperature. 
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It is interesting that the rate of increase in the central temperature in the third phase is 
smaller than that in the first one, whereas one might reasonably expect that the combined 
effect of both thermal diffusion and homologous compression would produce a more rapid 
increase in the central temperature during the third phase. The actual behavior may be 
explained by observing that the rate of energy loss via plasma neutrino emission increases 
with increasing temperature, thus leading to a reduced efficiency of heating of the central 
regions where temperatures are highest. 

In Fig. 9 are plotted, for every adopted accretion rate, temperature profiles at different 
times during the accretion process. For models with accretion rates larger than 10~^ 
Mq yr~^, the evolution of both temperature and density at the center is determined only 
by homologous compression. In fact, for such models, the thermal excess introduced by 
the accretion process into outer layers does not reach the center before runaway carbon 
burning (initiated at the point where the temperature is at a maximum) terminates the 
process. The model accreting matter at M = 2 x 10^^ Ai^ yr^^ defines a transition line: 
the physical conditions suitable for carbon burning are achieved first near the center (see 
panel d) in Fig. 9, M = 0.2 Mq), but quite soon also in outer layers where the maximum in 
temperature is located (M = 1.375 Mq). As time goes on, however, the rate of delivery of 
nuclear energy increases more quickly in the outer layers than in the deep interior and the 
thermonuclear runaway occurs first in the outer layers. 

EDITOR: PLACE FIGURE | HERE. 

In Fig. 10 are plotted the density profiles for the same models described in Fig. 9. 
The initial expansion phase is clearly shown by the decrease in the surface density whereas 
the action of the homologous compression is evidenced by the self-similar evolution of the 
density profiles. The same behavior is evident in the evolution of the temperature profiles 
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for the case M = 10~^ Mq yr~^; as time increases, the temperature level of the inner regions 
increases homologously while the continuous accretion process induces the formation of a 
peak in temperature in the outer layers. 



EDITOR: PLACE FIGURE 10 HERE. 



Models with 10""^ < M < IO^^Mq yr~^ {panels e and / in Fig. 9) achieve the conditions 
for carbon ignition in the inner zone when the maximum in temperature is still located in 
outer layers. This is also apparent from an inspection of Fig. 11, where we report profiles 
in the p — T plane for the models which have achieved ignition conditions. This behavior 
is due to the fact that the energy loss rate via plasma neutrinos, epiasma-i/ oc T^, is very 
large near the mass coordinate where the maximum in temperature is located. Hence, the 
gravitational energy delivered by the accreted matter causes a steepening of the temperature 
gradient and, paradoxically, the energy-loss sink due to ^piasma-i/ becomes deeper, reducing 
the amount of heat which can be transferred toward the center. Due to the effect of plasma 
neutrino losses, there is no carbon ignition for densities in the range 10^ — 10^ g/cm^. 



EDITOR: PLACE FIGURE 11 HERE. 



Models with accretion rates smaller than M < 10"'^ Mq yr~^ succeed in removing 
the energy excess deposited in outer layers because the gravitational energy is deposited 
more slowly than it diffuses inward. The surface energy reservoir is smaller and, hence, 
the energy losses via plasma neutrinos are smaller, so that the heat transfer inward is very 
efficient. In addition, the homologous compression induced by the accretion process is less 
rapid and the accreting star approaches the Chandrasekhar mass without carbon burning 
being ignited at a significant level anywhere within the structure. On the other hand, as 
the total mass approaches Mch , the rate of contraction becomes larger and the effect of 
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the induced compression, which is larger at the center, determines the physical conditions 
suitable for a central carbon ignition. During subsequent evolution under highly degenerate 
physical conditions, the energy delivered by carbon burning is stored locally, increasing the 
local temperature; the onset of convection (see panels / to i in Fig. 9) marks the onset 
of the thermonuclear runaway which will cause the disruption of the accreting star as a 
Type la Supernova. Although in these models the leading parameter in the evolution up 
to the carbon-ignition stage is the compression induced by strong contraction, the matter 
deposited on the surface plays a not completely negligible role in determining the ignition 
conditions. On decreasing the accretion rate, the density at which carbon ignition occurs 
becomes larger and, correspondingly, the temperature becomes smaller. 



4. The Effect of Rotation 

4.1. Physical Assumptions and Numerical Procedures 

To explore the effect of stellar rotation on the thermal response of CO WDs accreting 
matter we have included the lifting effect into the hydrostatic equilibrium equation 
according to the prescriptions in Dominguez et al. (1996): 

dP GM 



dr 



(15) 



where / is the ratio between the centrifugal and gravitational forces / = — = . In 

^9 9 

order to maintain the spherical symmetry in the treatment of the problem, the centrifugal 
force on a generic mass element dm is averaged over a sphere of radius r so that its value is 
provided by the following relation: 

Sr 47rr2 4 ^ ^ 

This treatment is equivalent to those by Kippenhahn & Thomas (1970) and Endal & 
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Sofia (1976, 1978) except that we neglect the corrective factor included by these authors 
in the expression for the radiative gradient (see Dominguez et al. 1996). Since WDs are 
chemically fairly homogeneous and are very compact objects, the efficiency of angular 
momentum redistribution may be large enough (Maeder & Meynet 2000) that the accreting 
WD can be treated as a rigid rotator, and we assume this to be the case. 

As discussed by Tassoul (1978), a rotating fluid undergoes instabihties which must be 
accounted for in a realistic description of the phenomenon. By adopting a rigid-rotator 
model, neither the Solberg nor the Goldreich-Schubert-Fricke instabilities arise. On the 
other hand, rotation could induce deformations of the equipotential surfaces in a such 
a way that meridional circulations on a large scale arise (von Zeipel paradox). In this 
case, deviations from spherical symmetry can be checked by computing a dimensionless 
quantity 7 defined as the ratio between rotational and gravitational energy (7 = Ej-ot/Eg). 
If 7 > 0.14, the equilibrium configuration is no longer spherical and meridional circulations 
cannot be neglected. In our computations, the 7 parameter is always smaller than 0.1, so 
meridional circulations have not been taken into account. Finally, if the centrifugal force 
acting on a mass element dm becomes larger than the gravitational force (/ > 1 in eq. 3), 
this element will experience a radial acceleration outward, achieving a velocity greater than 
the escape velocity (Roche instabihty) and the accretion process will be terminated. The 
angular velocity for which / = 1 is defined as the critical angular velocity uicrit, and its value 
is given by 



where r is the distance from the center to the mass element dm and g is the local value 
of the gravitational acceleration. This instabihty has been checked in detail during the 
evolution of the accreting models and will be addressed in the discussion of results. 

In the present computation, it is assumed that the accreted matter flows to the WD 




(17) 
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from a co-rotating accretion disk and that the deposited material has the specific angular 

momentum jo = uod?^ where d is the radius of the accreting WD and uo is the instantaneous 

angular velocity of the WD. The rate of change of the angular velocity uj of the WD model 

is then given by — = -M, where / is the moment of inertia of the model. 
at I 

The initial angular velocity, and, as a consequence, the initial total angular momentum 
of the accreting WD have been fixed parametrically by assuming that uj = a cUcrit.o, where 
i^crit,o is the initial value of Wcrit and a is a constant (< 1). In our case, cUcrit.o = 0.625 rad/s 
and, in the computations, we have adopted a = 0.1, 0.3, 0.5, 0.7, and 0.9. 



4.2. Numerical Results 

The inclusion of rotation in the hydrostatic equilibrium equation has a "lifting" effect, 
reducing the local value of pressure throughout the structure in such a way that the 
accreting star is both less dense and cooler than in the standard case (rotation neglected). 
While the density is only moderately decreased, the decrease in temperature throughout the 
accreting structure leads to a large enough increase in the thermal diffusion time scale to 
have a noticeable effect on the evolutionary behavior of the model. This is clearly shown in 
Fig. 12 where the initial value of Tdifr is plotted as a function of the initial angular velocity. 



EDITOR: PLACE FIGURE 13 HERE. 



On the other hand, the effect of compression induced by the accretion process is only 
slightly modified by the inclusion of rotation. Since Tdifr is far more sensitive to the inclusion 
of rotation than is Teh, an increase in the initial angular velocity of the WD always means a 
decrease in the value of the accretion rate below which central carbon ignition (as defined 
in §3.2) will occur. Evolution during the expansion phase at the beginning of the accretion 
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process is modified in a similar way. Since the initial expansion process is totally driven 
by the ratio Tdiff/Tch, for a fixed accretion rate (that is, for a given amount of gravitational 
energy delivered onto the WD per unit of time), the local energy excess in surface layers 
and, hence, the induced expansion become larger, the larger the initial angular velocity. 
This dependence is illustrated in Fig. 13, where the evolution of surface radius versus total 
mass is shown as a function of initial angular velocity. 



EDITOR: PLACE FIGURE ?? HERE. 



When rotation is not taken into account, the expansion comes to a halt due to a 
decrease in r^is in surface layers; gravitational energy is transferred efficiently to underlying 
zones from surface layers, which consequently begin to contract (see the discussion in 
§3). When rotation is taken into account, the Roche instability can occur, leading to an 
expansion phase during which tUcrit decreases in layers near the surface. In this case, the 
accreting structure may become gravitationally unbound after only a small amount of 
matter has been accreted. In Fig. 14 are shown the time evolution of several quantities for 
the model with M = IO^^Mq yr^^ and ujq =0.5 cc'crit.o- The quantities are the thermal 
diffusion time scale near the surface (panel a), the surface temperature (panel b), the radius 
(panel c). Panel d shows the evolution of the WD angular velocity uj (lower curve) and of 
the critical angular velocity in the outermost layers cUcrit.s- As the surface radius increases, 
i^crit.s decreases while, because of the deposition of angular momentum by accreted matter, 
UJ increases almost linearly with increasing WD mass. 



EDITOR: PLACE FIGURE 14 HERE. 



From these results, it follows that, for each value of M, there exists a threshold value 
ujth for ljq such that, if > Uth, the accreting structure becomes gravitationally unbound 
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due to the expansion during the first phase of the accretion process. In the M — ujq plane of 
Fig. 15, the shaded region gives the initial angular velocity and accretion rates for which 
the accreting WD becomes gravitationally unbound when its total mass approaches 1.4 Mq. 



EDITOR: PLACE FIGURE 15 HERE. 



Models which survive the expansion during the first phase of the accretion process 
recontract since the decrease in Xdiff in outer layers leads to an efficient transfer inward 
of the gravitational energy released. During this phase, the evolution of the accreting 
structure resembles the evolution of a non-rotating model, although, for a given value of 
the accretion rate, the whole structure is cooler and less dense for larger uoq . In Fig. 16, 
we show the same quantities as in Fig. 14 for a case in which the accretion rate is the same 
(M = IQ'^Mq yr~^), but in which the initial angular velocity is smaller (cuo/tUcrit.o = 0.1). 
In this case, the initial expansion is not large enough to significantly slow the decrease in 
TdifT, and both uj and ujcv\t,s increase steadily with increasing WD mass as the energy excess 
deposited in outer layers is efficiently transferred toward the center, as is clearly depicted 
in Fig. 17 where temperature profiles for various models during the accretion process are 
shown. 



EDITOR: PLACE FIGURE 16 HERE. 



As detailed in Fig. 16, the accreting WD spins up as its radius increases, with both 
UJ and cUcrit.s increasing steadily; as a result, the temperature level throughout the WD 
increases homologously (see also Fig. 17). When the WD mass reaches ~ 1.3 Mq, the 
centrifugal force begins to compete with the gravitational force; thereafter, the increase 
in the total angular momentum of the structure causes the surface layers of the WD to 
expand and the compression induced by the accreted matter is no longer able to balance 
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the radiative losses. As a consequence, the entire structure cools down (see the last profile 
in Fig. 17)! Notice further that the expansion produces a decrease in Ucrit to the extent 
that the accreting structure becomes gravitationally unbound; the accretion process comes 
to a halt. 



EDITOR: PLACE FIGURE 17 HERE. 



It thus appears that the occurrence of the Roche instability prevents the accreting 
models from achieving the physical conditions necessary for the ignition of carbon; the 
accreting structures are not only too cool but, in addition, the centrifugal force due to 
rotation prevents the contraction in interior regions which produce heating by the conversion 
of gravitational potential energy. 



4.3. The Dependence of the Final Outcome on jd 

The scenario summarized in Fig. 15 depends on the amount of angular momentum 
deposited by the accreted matter: in fact, the increase in the angular velocity, particularly 
during the early phase of the accretion process, is totally controlled by the deposition 
of angular momentum. In order to make this point clearer, we have re-computed the 
previously discussed models assuming that the accreted matter flows to the WD from an 
accretion disk at rest (jd = 0). In Fig. 18, we report the evolution of the surface radius 
for models accreting CO-rich matter at the rate M = 10~^ Mq yr~^ with different initial 
angular velocities, assuming that the accreted matter deposits (dashed lines) or does not 
deposit (solid lines) angular momentum. It is clear that, if angular momentum is not added, 
the accreting structure with initial angular velocity smaller than 0.9 cUcrit.o survives the 
initial expansion and thereafter contracts as a consequence of the increase in the total mass. 
In addition, if jd = 0, for uo in the range 0.3 ~r 0.9 co'crit,05 the accreting structure ignites 
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carbon burning at the center since the total mass is larger than 1.4 Mq. It is also clear that, 
the larger the initial value of the angular velocity, the larger is the final total mass of the 
accreting WD. This occurrence is consistent with the known fact that the Chandrasekhar 
mass for rotating structure is larger than that for non rotating models, depending strongly 
on the angular velocity (see Anand 1968, Ostriker and Bodenheimer 1968). 



EDITOR: PLACE FIGURE M HERE. 



By assuming that the accreted matter does not deposit angular momentum, the 
morphology of the M -ujq plane changes from that shown in Fig. 15 to that shown in Fig. 
19. The changes are worth emphasizing: 

• for large values of both ujq and M (the region labeled "Mj < 1 Mq"), the accreting 
models experience the Roche instability when only a small amount of matter has been 
accreted (AMacc < 0.1 M©). This zone corresponds to the non-shaded region in Fig. 
15, so that its lower boundary represents the cuth previously defined; 

• for large values of c^o and small values of M (the region labeled "^n < ^v^), the 
accreting structures never achieve the physical conditions necessary for carbon 
ignition, although they do not become gravitationally unbound. The computation of 
these models has been halted because the adopted input physics, in particular the 
equation of state, is inadequate. In any case, the final mass of the accreting structure 
is larger than 1.4 M©; 

• for small values of both ujq and M (the region labeled "Central Ignition"), the 
accreting structures ignite carbon burning at the center when the total mass 
approaches (and in some case exceeds) the non-rotating Chandrasekhar mass; 



- 25 - 



• the zone labeled "Central and Off-Center C-ignition" represents structures which 
ignite carbon burning both at the center and in the outer layers. This zone marks the 
transition between models which become gravitationally unbound and models which 
succeed in igniting carbon burning; 

• models in the region above the dashed line and below the present Uth line become 
gravitationally unbound when the total mass of the accreting structure approaches 
the non-rotating Chandrasekhar mass. This region corresponds to the shaded zone in 
Fig. 15. 



EDITOR: PLACE FIGURE 19 HERE. 



It is evident that, on including the lifting effect of rotation on the evolution of a CO 
white dwarf accreting CO-rich matter, the zone in the Mwd - M plane suitable for a 
central carbon ignition becomes narrower. If matter is accreted at rest, carbon burning 
is ignited at the center for the same values of the accretion rate as in the standard case 
(M < IO^^Mq yr^^), but only for small values of the initial angular velocity (see Fig. 19), 
whereas, if the matter is accreted from a co-rotating disk, the WD never ignites carbon 
burning. In this latter case, the increase of the angular velocity due to the deposition of 
angular momentum by the accreted matter induces the expansion of the accreting structure 
which becomes gravitationally unbound. 



4.4. The Braking Phase 

Interestingly enough, rotating models with small values of both the initial uj (namely 
uq < 0.3 tUcrit.o) and accretion rate (namely M < 10^^ Mq jt^^ ) undergo the Roche 
instability when their total mass is of the order of the non-rotating Chandrasekhar mass. 
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In this case, the ratio of rotational and gravitational energies (7 = Ej-o^/Eg) approaches 
the critical value ~ 0.14 and, as a consequence, the accreting structure has to deform into 
an ellipsoidal shape (see §4.1); due to the induced deformation, the mass distribution of 
the white dwarf becomes anisotropic in such a way that gravitational wave radiation is 
emitted (see Shapiro & Teukolski 1983, 1990). The emitted waves carry away rotational 
energy and angular momentum and the total angular momentum of the braking structure 
evolves according to an exponential law J — Joe~*/^, where Jq is the value of the angular 
momentum at the beginning of the slowing down phase and r is the time scale on which the 
process occurs; the latter quantity is determined by the efficiency of braking and, therefore, 
by the angular velocity of the structure when 7 approaches the critical value. Since the 
intensity of the emitted GWR depends on o;^, we assume that the braking process becomes 
the driving mechanism for the evolution of the accreting white dwarf when the angular 
velocity is of the order of 1 rad s~^. 

It is worth noticing that, on decreasing the total angular momentum, the effective 
gravity acting on a given mass element increases and, as a consequence, the whole structure 
experiences a contraction which produces the heating-up of the entire star. In addition, 
the contraction induced by the braking produces an increase in the angular velocity of the 
white dwarf so that the 7 parameter actually remains of the order of the critical value. This 
means that the white dwarf continues to emit GWR copiously as it contracts. 

These considerations suggest that the braking process acts as a compression mechanism 
and that the r parameter in the exponential braking law can be interpreted as a 
compressional heating time scale; we can therefore argue that the final outcome of the 
braking phase depends on the ratio between r and Tdifr, the thermal diffusion time scale of 
the considered structure. 

In order to explore this conjecture, we have constructed an evolutionary sequence 
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which begins with a model at the end of the previously explored sequence with initial 
angular velocity Uq = O.lcUcnt.o = 0.067 rad/s and accretion rate M = 10^^ Mq yr^^ . In 
this last model, 7 is of the order of the critical value, the total mass of the white dwarf is 
~ 1.42 Mq , the angular velocity is ~ 1.1 rad/s, and the value of r^is is ~ 2 x 10'^ yr; at 
this time the structure is isothermal, as shown in Fig. 17. 

In order to explore the effect of GWR braking on further evolution, we choose two 
values of r: r = 10'' yr, which corresponds to a low braking efficiency, and r = 10^ yr, 
which corresponds to a high braking efficiency. As discussed in Shapiro & Teukolski (1983) 
and in Ma (2002), the GWR-driven braking process is not important on time scales greater 
than 10^ yr since, then, the slowing down is dominated by viscosity. 

In the first case (r = 10'' yr), the homologous compression induced by GWR braking 
produces an increase in temperature everywhere in the structure, as shown in Fig. 20, 
where temperature profiles for selected models are presented. When the physical conditions 
suitable for carbon burning ignition are achieved at the center, the released nuclear energy 
is stored locally, due to the high electron degeneracy; a convective core forms and grows 
rapidly until convection reaches the surface. The thermonuclear runaway will ultimately 
proceed on a dynamical time scale, causing the disruption of the white dwarf as a type I 
supernova. 



EDITOR: PLACE FIGURE 20 HERE. 



In the second case (r = 10^ yr), the GWR- induced compression also produces a 
homologous increase in the temperature level throughout the structure. But, superimposed 
on this increase is an additional increase in layers near the surface, where a local temperature 
maximum forms (see Fig. 21). This secondary maximum occurs in part because matter near 
the surface, being less electron degenerate, has a larger specific heat. Additionally, however. 
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the evolution of the model occurs on a time scale smaller than the thermal diffusion time 
scale in surface layers, so the thermal energy released into layers near the surface remains 
localized in these layers. In any case, the evolution of the model is totally dominated by 
the homologous compression due to GWR-induced braking, and, just as in the first case, 
the thermonuclear runaway which has begun at the center will evolve into an explosion of 
supernova proportions. 

EDITOR: PLACE FIGURE |2T| HERE. 

Thus, the final outcome in both cases is a SN la event, independent of the assumed 
GWR-braking efficiency. The only effect of a variation in the braking efficiency over a 
substantial range is a modest variation in the physical characteristics of matter when carbon 
burning first develops into a thermonuclear runaway at the center. This is made clear in 
Fig. 22, where we show the evolution of central characteristics in the p — T plane for the 
two discussed models. 



EDITOR: PLACE FIGURE 22 HERE. 



5. Final Remarks 

In this work, we have explored the lifting effect of rotation on the evolution of a CO 
WD accreting CO-rich matter from the thick disk that is initially formed during the merger 
of two close binary CO WDs. We find that rotation induces a cooling of the accreting 
structure to the extent that the final outcome of the accretion process can be significantly 
modified relative to the case when rotation is neglected. For a fixed value of M and an 
initial angular velocity of the WD larger than a critical value, the accreted matter produces 
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a large enough thermal energy excess in outer layers that the star itself expands to a larger 
radius. The induced expansion leads to a reduction of the surface gravity to the extent that 
the accreting structure becomes gravitationally unbound (Roche instability) and no more 
matter can be added. We have shown that this occurrence depends on both the accretion 
rate and the initial angular velocity. For large values of both M and ojq , the accreting 
WD becomes gravitationally unbound when only a small amount of mass has been added 
(AMacc < 0.1 Mq), whereas, for small values of both M and uq , the WD attains a total 
mass close to 1.4 Mq before the Roche instabihty sets in. We have also shown that behavior 
of this sort depends critically on the angular momentum carried by the accreted matter. 
For a fixed accretion rate, the larger the specific angular momentum of the accreted matter, 
the earher does the Roche instability occur (the WD spins up more rapidly) . 

The behavior we describe is a consequence of the fact that the thermal diffusion time 
scale is very sensitive to changes in the assumed rotation rate (increasing with increasing 
rotation rate) whereas the compressional heating time scale is almost independent of such 
changes; the rate of transfer toward inner regions of the energy excess produced in outer 
regions by accreted matter becomes smaller and smaller with respect to the standard case 
(without rotation) as the rotation rate is increased. As a consequence, for a given accretion 
rate and WD mass, the density and temperature at the center are smaller, the larger the 
initial angular velocity. In addition, we find that, if the specific angular momentum of the 
accreted matter is equal to that in the outer layers of the WD, all the computed models 
become gravitationally unbound well before the physical conditions suitable for carbon 
burning ignition are attained anywhere in the WD structure. 

For small values of both initial angular velocity and accretion rate (namely, 

ujq < 0.3cJcrit.o =^ 0-18 rad/s and M < IO^^Mq yr^-*^ ), the accreting model becomes a very 
fast rotator (a;f ~ 1 rad/s) when the total mass is of the order of 1.4 Mq and 7 (= Ej-ot/Eg) 
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becomes of the order of the critical value (7^ ~ 0.1). Then, the accreting WD adopts an 
elliptical shape and begins to emit gravitational waves which carry away both angular 
momentum and rotational energy. The important point is that the braking induced by 
GWR produces a contraction of the whole structure and, as a consequence, the WD spins 
up in such a way that the 7 parameter continues to be of the order of the critical value. 

In this way, the evolution of the structure is totally driven by the homologous 
compression induced by the GWR-braking process. According to our results, independently 
of the efficiency of braking, which means independently of the time scale on which braking 
occurs, the physical conditions suitable for an explosive carbon-burning ignition are attained 
at the center of the WD, so that a SN la event occurs. 

For models with large cuq and M , the Roche instability occurs when the total mass of 
the WD is smaller than 1.4 Mq , the exact value depending on the initial values of both the 
accretion rate and the angular velocity. But this result has been obtained by assuming that 
the accretion rate is strictly constant; in the real world, the WD will continue to accrete, 
but at a smaller rate. As the Roche instability is approached, the accretion rate drops, the 
angular momentum of the WD increases less rapidly, and the rate of producing an energy 
excess in the surface layers decreases. Thus, a self-regulatory mechanism prevents the 
gravitational instability from actually taking place. Rotation acts as a tuning mechanism 
for the accretion rate in such a way that the accreting structure can continue to accept 
mass even in the case of a very large initial accretion rate. We address the details of this 
topic in a forthcoming paper. For now, it is sufficing to observe that, although the accretion 
rate decreases due to the lifting effect of rotation, the total angular momentum and, as a 
consequence, the angular velocity of the WD must increase; the WD will become a fast 
rotator. Gravitational wave radiation becomes competitive with thermal diffusion when 
the angular velocity of the WD attains a value of the order of 1 rad s~^, and, thereafter. 
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GWR-induced braking leads to the physical conditions suitable for carbon-burning ignition 
in the center and, hence, to a SN la event. 
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Table 1. Results of experiments with standard conditions on a O.SMq CO WD accreting 
at different accretion rates. M is tfie accretion rate (in Mq yr""^), Mf (in Mq) is tfie total 
mass at carbon ignition, Mig is the mass coordinate where carbon burning is ignited (in 
Mq), Tig (in K) and pig (in g cm^^) are, respectively, the temperature and the density at 

the ignition point. 



M 


Mf 


Mig 


log(rig) 


log(Pig) 


10-5 


1.037 


0.899 


8.785 


5.977 


5 X 10-6 


1.193 


1.134 


8.773 


6.242 


3 X 10-6 


1.315 


1.297 


8.772 


6.390 


2 X 10-6 


1.383 


0.000=^) 


8.295 


9.565 


10-6 


1.380 


0.098 


8.381 


9.371 


5 X 10-7 


1.375 


0.000 


8.348 


9.392 


10-7 


1.380 


0.000 


8.287 


9.505 


5 X 10-^ 


1.383 


0.000 


8.237 


9.576 


io-« 


1.395 


0.000 


8.101 


9.705 



During the subsequent evolution of this model, the mass coordinate where the nuclear 
energy production is maximum jumps to the point M = 1.375 M© just below the surface. 
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Fig. 1. — Chemical profiles of ^^C and ^^O in a 0.8 Mq model at the beginning of the 
cooling sequence as obtained by evolving a pure He star (panel a) and as it becomes after 
the re-homogenization made according to Salaris et al. (1997) [panel b) (see text). 

Fig. 2. — Time evolution of the thermal diffusion time scale (solid lines) and of the 
compressional heating time scale (dashed lines) for all of the computed models as listed 
in Table 1. 

Fig. 3. — Time evolution of the central temperature (sohd line) , the thermal diffusion time 
scale (dashed line), the mass coordinate where the maximum temperature is located (long 
dashed line) for the model with M = 10"^ Mq yr~^ (see text). All the quantities have been 
normalized to their maximum value. 

Fig. 4. — Evolution of the temperature and of ratios of several time scales in the interior 
of a model accreting at the rate 10~^ Mq yr~^. Here, Tdiff,ioc is a local diffusion time scale, 
7'comp is a local compressional time scale, and Tevo is the time elapsed from the beginning of 
the accretion experiment. 

Fig. 5.— Same as Fig. 4, but for M = 10~^ Mq yr^^ 

Fig. 6. — Evolution in the HR diagram of all of the computed models hsted in Table 1 (see 
text). Attached to each line is the corresponding value of M (in Mq yr~^ ). 

Fig. 7. — Time evolution of the surface radius for all of the models listed in Table 1 (see 
text). Attached to each line the corresponding value of M (in Mq yr~^ ). 

Fig. 8. — Evolution in the p — T plane of the center of the 0.8 Mq CO WD accreting matter 
at various M, as labeled inside the figure (in Mq yr~^ ) (see Table 1). The long-dashed curve 
labeled " e„ — Scc " is the ignition curve. 
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Fig. 9. — Temperature profiles for several selected structures during the accretion process 
for different values of the accretion rate, as labeled inside each panel. 

Fig. 10. — Density profiles for the same models as in Fig. 9. 

Fig. 11. — Profiles in the p — T plane of models which have attained the physical conditions 
suitable for carbon ignition as a consequence of mass accretion at different M, as labeled 
inside the figure (in units of Mq yr~^). 

Fig. 12. — Initial value of the thermal diffusion time scale Tdiff (in yr) as a function of the 
initial angular velocity (in units of a;crit,o) ■ 

Fig. 13. — Time evolution of the surface radius for models rotating with different initial 
angular velocities as labeled inside the figure (in units of Ci;crit,o)- Each panel refers to a 
different value of the accretion rate as labeled. For comparison the evolution of the standard 
model without rotation (dashed line) is also shown. 

Fig. 14. — Evolution of the diffusion time scale in the outer zone {panel a), the surface 
temperature {panel b), and the total radius {panel c) of the model with M = 10~^ Mq yr~^ 
and cuq — 0.5a;crit,o- In panel d, the evolution of the white dwarf angular velocity (lower 
curve) is compared with that of the critical angular velocity at the surface (upper curve) . 

Fig. 15. — The region in the parameter space M — loq (shaded zone) for which the rotating 
models accreting CO-rich matter become gravitationally unbound when the total mass of 
the WD approaches 1.4 Mq (see text). 

Fig. 16. — The same as in Fig. 14, but for the case M = 10~^ Mq yr~^ and ujo = 0.1a;crit,o- 

Fig. 17. — Temperature profiles of selected structures during the evolution of the rotating 
model with luq = 0.1a;crit,o and accreting matter at M = 10"*^ Mq yr~^. The labels attached 
to each line provide the time elapsed from the beginning of the accretion process. 
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Fig. 18. — Evolution of the surface radius for models accreting CO-rich matter at M = 
1O~^M0 yr~^ with different values of the initial angular velocity computed by assuming that 
the accreted matter does (dashed lines) or does not deposit (solid lines) angular momentum 
(see text) . 

Fig. 19. — The same as in Fig. 15, but for rotating models which accrete CO-rich matter at 
rest (see text) . 

Fig. 20. — Temperature profiles for some selected structures during the braking stage of a 
rotating ~ 1.42 Mq CO WD. The braking time scale has been set equal to 10^ yr (see text). 

Fig. 21. — Temperature profiles for some selected structures during the braking stage of a 
rotating ~ 1.42 Mq CO WD. The braking time scale has been set equal to 10^ yr (see text). 

Fig. 22. — Evolution in the p — T plane of the center of models accreting CO-rich matter at 
M = 10^^ Mq yr^^. The dashed line describes the evolution of the standard model (a; = 0), 
whereas the solid lines describe the evolution of two rotating models with cuq — 0.1a;crit,o and 
two different GWR- braking times, one with r = 10^ yr, the other with r = 10^ yr (see text). 



7 



6 



7 



6 



7 



6 



5 



1 1 1 1 1 1 1 1 1 1 1 1 1 

L a) : 

" M=10-5 

pTTT| 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 

- d) : 

L M=2xl0-8J 


1 1 1 1 1 1 1 1 1 1 1 1 1 

L g) M=io-^ J 


- b) : 

M=5xl0-6_ 


1 1 1 1 1 1 1 1 1 
- e) J 

L M=10-6 J 


1 1 1 1 1 1 1 1 1 

L h) M=5xl0-8J 

1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 

- c) : 

L M = 3xl0-6J 

1 1 1 1 1 I 1 I 1 1 1 i 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 
- f) - 

L M = 5xlO-^J 
1 I 1 1 1 1 1 1 1 1 1 1 1 


/ i 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.8 1 1.2 1.4 0.8 1 1.2 1.4 0.8 1 1.2 1.4 

M/Mq M/Mq M/Mq 



4 



2 







-2 



1 I 1 1 1 1 1 1 1 1 1 1 

10-4^ 


1 1 1 1 1 

— 


^ VsxlO-s 








- 1 n-6 
<,X iU X 




10-6 X 


- 


5x10-' \^ 




10-' 




5x10-3 \ 


— 


10-8*~"~~~~--^.^^^\ 




1 I 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 



5.5 5 4.5 4 



log(Tj 




log(M) 




7 8 9 10 



